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Abstract 

We describe a 8085 microprocessor interface developed to make reliable time period 
measurements. The time period of each oscillation of a simple pendulum was measured 
using this interface. The variation of the time period with increasing oscillation was 
studied for the simple harmonic motion (SHM) and for large angle initial displacements 
(non-SHM). The results underlines the importance of the precautions which the students 
are asked to take while performing the pendulum experiment. 
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1 Introduction 

Qh| The simple pendulum is a very trivial experiment that physics students do in higher secondary. 
Yet students sometimes fail to appreciate why "the initial angular displacement of the pen- 
dulum must be small" . Interesting letter to the Editor in Physics Education (UK) point that 
\q . the pendulum's time period increases by only 1% for pendulum's oscillating through 30°. A 
1% increase means the time period is 10msec more for than a pendulum undergoing SHM 
whose time period is lsec, thus students do question the very need of the precaution that 
the pendulum should only be given small angle initial displacements. Variations as small as 
10msec are very difficult to measure using stop-watches. Since computers (or microprocessors) 
can in principle make measurements in micro-seconds, we were tempted to study the simple 
O ! pendulum using a micro-processor. 

2 The interface 
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The microprocessor is essentially made up of digital devices, which communicate among itself 
in the language of binary, i.e. in ones (1) and zeros (0). That is voltage signals of preferred 
levels. To find the time interval of an oscillating pendulum, we keep an arrangement of laser 
source and light dependent resistor (LDR) such that the pendulum's bob cuts the light path 
during its oscillation. As the bob cuts the path, the light is momentarily blocked. This 
produces a change in current generated in the LDR. For the microprocessor to communicate 
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and understand this change (analog) in current, it has to be converted to TTL compatible 
digital voltage. The conversion and subsequent wave shaping is done using the circuitry shown 
in figfllj. 
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Figure 1: Circuit used to wave-shape and interface instant of pendulum's cross-over in front 
of the light dependent resistance (LDR). 

The value of resistance of the LDR, as the name suggests, depends on whether light is falling 
on it or not. The resistance of the LDR is usually inversely proportional to the intensity of 
light falling on it pQ. In our case when the laser light was falling on the LDR, it's resistance 
was 15KQ, while on switching off the laser light it rose to 150Kfi. The voltage drop across 
the resistance which forms a voltage divider with the LDR is 5v when the laser light is 

'ON' (bright phase) and 0.75v when LDR is not exposed to laser (dark phase). This voltage 
(V i) could have been directly fed into the microprocessor, however, the intensity of the laser 
light strongly depends on the current supplied by it's batteries. With time since the current 
is likely to fall, the voltage V i, would change as the experiment is being conducted. The 
difference amplifier, amplifies the difference between the voltages V i and Vi. By selecting a 
proper Vi, using the pot, the output V o2 varies from positive level for dark phase and negative 
level for the bright phase. This inversion is bought about by the inverting amplifier (opamp). 
The second opamp inverting action brings V Q 3 in phase with V i. This opamp is essentially a 
Schmitt's trigger, which hard drives the output to +12v and -12v. The output of the second 
opamp (V 3) varies from -12v to +12v in accordance with the motion of the pendulum. Since, 
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the microprocessor can only understand zero or high state (~ +5v), a 4.7v zener diode is used 
to protect the the microprocessor (from +12v) by converting + 12v level to +4.7v and to force 
-12v level to zero. The 1KQ resistance kept between the opamp and zener diode is to control 
the current flowing into the zener diode. 
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Figure 2: Plot of pendulum length vs time period squared (T 2 ). The data point fall on a 
straight line with co-relation factor as good as 0.9995. 

The realibilty of our microprocessor program (listed in the Appendix) was checked by 
finding the frequency of known square waves fed from audio function generators. The program 
is essentially a counter program which counts the time interval taken between two positive 
edges of a train of square waves. While counting, the program loops between instruction 
addresses C024 H and C02D H for the square wave's high state and between C02E H and C037h 
for the input square wave's ground state. The latter part was to over come the requirement 
of IC555 monostable trigger hardware in our circuitry. For a count of 'N', the time period is 
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given as 



40iV - 3 

= — 

/AON - 3 
V 3 



fis (1) 



where is the frequency of the microprocessor clock in MHz. The formula is obtained by 
calculating the 'T' states or time taken for the microprocessor to execute each instruction. It 
should be noted, the microprocessor's program takes 17/isecs to identify an edge and write the 
count in a memory location. This 17/xsecs is a systematic error that would be present in the 
value of time period measured. From this exercise we realized the quartz crystal used as a 
clock for the microprocessor kit was 6.2MHz (Books j2] state it is 6MHz, however, it depends 
on local manufacturers). Thus, results of all time measurements (fig 2, 6, 8 and 11) have to 
be multiplied by 0.968(=6.0/6.2) to accommodate this correction. 

While the realibility of the software was established, to make sure the LDR (our transducer), 
as also our microprocessor interface was reliable, we measured the time period of oscillation 
for various pendulum lengths. Figure (J2J) shows a perfect linearity between the pendulum's 
length and it's time period squared (i.e. T 2 ). This is expected and is in accordance to well 
established theory that we shall discuss below. It should be noted that on an average the 
response time of a LDR is in small milli-secs. Thus, the ability to resolve and measure any 
changes in time period with increasing oscillations would be in milli-secs. Also, the accuracy in 
terms of absolute value of 'g' calculated from experimental data would depend on how accurate 
the quartz crystals frequency has been reported. 



3 Simple Harmonic Motion: small initial displacement 

Before proceeding to discuss the results of our experiment it would be worthwhile to recapulate 
about the pendulum and under what conditions does it's motion reduce to a simple harmonic 
motion. A pendulum is easily set up by suspending a point mass. Physically, this is achieved 
by suspending a bob which has an appreciable mass but whose radius is small as compared to 
the length of the string used to suspend the bob in consideration. The pendulum is set into 
to and fro motion by displacing it from it's mean position. The forces acting on the displaced 
pendulum is shown in fig®- The restoring force is given as 

F = — mgsinO 

where 'm '& and 'g ' are the mass of the bob, it's angular displacement from the mean position, 
the acceleration due to gravity respectively. The above leads to the equation of motion 

m L—— = — mqsvn<o 
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Figure 3: Simple Pendulum of length L. Also shown is the plot of error (in %) on assuming 
sin0 = G increases with increasing G. 



or 



On considering the initial angular displacement i.e. G to be small, sinQ of eqn(j2J) reduces to 
G and substituting u 2 = g/L, we have 

IF = ~ U 9 (3) 

This second order differential equation describes the motion of the simple harmonic motion 
(SHM), whose analytical solution is easily derivable and given as 

Q(t) = A sin{ut) + B cos{ut) 

where A and B are constants. We can get the values of the constants by choosing suitable initial 
conditions. The time period of oscillation can be obtained from the relationship {u = |?) 
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The above equation shows the proportionality between and the pendulum's length. To 
confirm the reliability of our time measuring device (interface and software etc), we confirmed 
this relationship, see fig(j2J). This relationship holds true for small angle displacements. Hence, 
the data for fig(J2J) were collected for various lengths of the pendulum with initial angular 
displacement being 5°, an universally accepted small angular displacement. 

Students identify eqn(jU easily, since it is used by them to estimate the acceleration due to 
gravity. Also of interest is the fact that the above expression implies that the time taken to 
complete one oscillation is independent of the angular displacement (6), abide subjected to the 
condition sin© ~ G. It's here that the argument starts as to what would be the appropriate 
initial displacement that a experimentalist should give to attain the simple harmonic motion? 
As (in radians) increases, the disparity between itself and it's sine (sinG) increases. This 
fact is seen in fig©, where the increase in disparity is shown in terms of error ( Q-gnQ ? j n 
w.r.t. G. 

As can be seen, the error is below 10% for angles less than 45°. Would this limit be 
acceptable? Before answering this question, as to understand the boundary between SHM 
and non-SHM, we proceed to understand the modifications introduced in eqn(@J), when the 
pendulum is set into motion with large angle displacements (non-SHM). 

4 Pendulum with large initial displacement 

The time period of oscillation of a pendulum oscillating with large angles can be found by 
solving eqn(J2J), i.e. 



However, that is easier said then done. Infact discussions on large amplitude oscillations are 
rarely carried out because there are no analytical solutions for the above differential equations. 
Infact, the solution is expressed interms of elliptical integrals [211] 



d 2 e 



— oj 2 sinQ 



(5) 



dt 2 





Hence, eqnfjHJ) is either numerically solved or various approximations are used. Of these ap- 
proximations, the most famous was given by Bernoulli in 1749 [4J 



where T Q is the time period had the SHM condition been satisfied and is given by eqn(flj) 
and m is the maximum angular displacement given to the pendulum. Eqn(jZJ) would suggest 
that take whatever initial displacement you want while doing the experiment to determine the 
acceleration due to gravity, all you have to do is to include the correction (0^/16) in the time 
period expression (eqnEJ). 
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Figure 4: Plot of acceleration due to gravity with initial displacement. The solid line is essen- 
tially calculated from 'T' of eqntJZj) while the constant line is 'g' evaluated using T Q = 1+e T ^ 16 
(i.e. time period after correcting for large angle displacement). 

While the pendulum and it's time period of oscillation in itself is interesting, it is usually 
used to evaluate the acceleration due to gravity. Teachers insist that students do the experiment 
with small angular displacements. Students do this without appreciating "why" and as to 
"what is a small angle". Let us consider what kind of variation is expected theoretically if this 
precaution is not adhered to. Figure(jH) shows the plot of the variation of acceleration due to 
gravity with the initial displacement done by numerically solving eq(jSJ) (Ej- The values of 'g' 
represented by the solid line is essentially calculated from 'T' of eqn(0) while the constant line is 
'g' evaluated using T Q = 1+e T j 1Q (i.e. time period after correcting for large angle displacement). 
As the figure shows, with increasing angular displacement, the error in evaluated 'g' grows. 

EqnfJZJ) suggests only a trivial consideration of including a correction factor is required if 
the small angle precaution is not followed, the question then arises "why fuss over small initial 
displacements?." Also, it is evident from fig(@J) that the error in 'g' would be below 10% for 
initial displacements below 45° which is quite a large angle. This might be well within the 
limits of experimental error, induced by your measuring devices like scale and stop-watch. It 
should be noted that Nelson and Olsson j^j have determined 'g' with an accuracy of 10~ 4 using 
a simple pendulum by including as many as 8-9 correction terms. Thus, the importance of 
maintaining small initial displacement while performing the experiment is still not convincing. 
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In the next section we report the measurements made by our microprocessor interface and try 
to address the questions we have asked above. 



Figure 5: Important instances for measuring one complete oscillation, with reference to our 
experiment. 

A bob of radius 2.5cm was suspended using a cotton thread of length one meter (the length 
of the pendulum thus is 1.025m). As the pendulum cut the lasers path to the LDR, an electric 
pulse is generated. From the point of the onset of the positive edge, the microprocessor counts 
the time elapsed till the bob cuts the light path again. In one complete oscillation, the bob 
cuts the light path thrice, say at instances t D , ti and t 2 (see fig EI)- The time period is given 
as t 2 — t D . The program is however designed to store ti — t D and t 2 — ti- This was done to 
make sure there is no error induced due to the inability to pin point the mean position of the 
pendulum. Data were collected for the pendulum oscillating with the initial displacements of 
5°, 10°, 15°, 20°, 25° and 30°. 

Fig® shows the variation of the pendulum's time period with oscillations. While for small 
initial angles, there is no or slight variation in time period, for large initial displacements, 
namely 20°, 25° and 30°, the fall in time period with increasing oscillations is pronounced. 
The fall is better appreciated by plotting the normalized data (T/Ti). Deviation from the 
smooth variation in time period (scattering of data points) is seen for large angles. This 
primarily is due to the pendulum's support not being perfectly stationary. For this reason we 
restrict our report to maximum angular displacement of 30°. 

The fall in time period seems exponential. Since the experiment is not ideal, one can expect 
damping to play a role (the decrease in amplitude was visible to the naked eye with increasing 



5 Results & Discussion 
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Figure 6: Variation of the pendulum's time period with oscillations as measured by the mi- 
croprocessor interface. Along side is the normalized (T/Ti) time periods variation with oscil- 
lations. 



number of oscillations). Infact, damping is expected to attenuate the amplitude of oscillation 
exponentially with time. Thus, eqn(0) can be written as 

T = T Q M + (8) 

Or, can be written as 

T = a (l + be~ m ) (9) 

Our objective would be to fit the above the equation to the experimental data of fig©- Table 
I lists the coefficient of eqntjUJ) obtained by curve fitting. 

Table I: List of the coefficients obtained by fitting eqn( |9"j) to the experimental data of fig©. 
The last column lists the co-relation of fit with respective data points. 



S.No. 


e 


a 


b 


P 


T(/3 = 0) 


r 




(degrees) 


(sec) 


(radian) 2 


xl0~ 2 sec _1 


(sec) 




1. 


5 


2.00441 


0.000333 


0.34751 


2.00508 


0.8138 


2. 


10 


2.00676 


0.001552 


0.31740 


2.00988 


0.9911 


3. 


15 


2.01014 


0.002072 


0.38300 


2.01430 


0.9968 


4. 


20 


2.00181 


0.006540 


0.32090 


2.01490 


0.9996 


5. 


25 


2.00519 


0.009755 


0.41942 


2.01038 


0.9921 


6. 


30 


2.00961 


0.011860 


0.50847 


2.03344 


0.9976 
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Figure 7: Plot of T(/3 = 0) (or a+ab) vs angle. The data points fit well to Bernoulli's approx 
(eqnEj) with co-relation factor as good as 0.985. 



Instead of using coefficient 'a' as a variable for curve fitting, it should be taken as a constant 
(T , the SHM time period) defined by eqn(jH). That is, for all six values of 9 (Table I), the 
value of 'a' should work out to be the same. However, this proved to be a difficult exercise 
where we were not able to achieve good co-relation between the experimental data points 
and the fitted curve. To overcome this, after obtaining the generalized coefficients (a, b) 
we calculated T(/3 = 0) (=a+ab, listed in Table I), i.e. the time period of the pendulum's 
oscillation through large angles without any damping (the variation with angle of oscillation is 
given by Bernoulli's approximation eqnEJ). Fig( | 7 jl shows the variation of T(j3 = 0) with angle. 
The solid line shows the curve fit (eqn[7J). The T Q determined from our experimental data 
works out to be 2.00545sec. One can use this value in eqn(jU) and determine the value of 'g', 
the acceleration due to gravity. We get the acceleration due to gravity as 9.963m/sec 2 . This 
value is on the higher side. Even after accounting for the influences corrections terms discussed 
in reference [E], the value of 'g' would be on the higher side. The error, might be due to the 
lack of precise knowledge of the micro-processor's clock frequency. Using eqQ and eq(jU, we 
get the maximum possible error as 



For an error of +0.1 x 10 6 MHz is the knowledge of the micro-processor's clock frequency, the 
maximum error in 'g' would be +0.65m/s 2 . Thus the error, as stated, is due to the imprecise 
knowledge of the micro-processor's clock frequency. 

We now investigate the remaining coefficients 'b' and (3. The coefficient 'b' is proportional 
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Figure 8: Plot showing variation of coefficient 'b' with increasing initial angular displacement. 
The solid line shows that the experimental data points fall on a parabola. Second plot shows 
how the damping coefficient of the pendulum varies with initial angular displacement. 



to G 2 (compare eqn |H] and EJ) . This is evident from fig© which shows the data points to fall 
nearly perfectly (co-relation factor is 0.984) on a parabola. The proportionality constant by 
eqn(jHJ) should be 1/16 (if 6 is in degrees), or 1.9 x 10~ 5 rad 2 . Our result gives the proportion- 
ality constant as 1.4 x 10~ 5 rad 2 (or 1/21). EqnfjHJ) is only an approximation hence, we can 
confidently say that our data follows the solution given by Bernoulli (eqnEJ). 

The second graph of FigfjHJ) shows the plot between the damping coefficient and the angular 
displacement. The damping factor is nearly constant for small angular displacements which 
shows a rapid linear increase for angular displacements above 20°. As can be seen the variation 
is similar to that of the IV characteristics of a diode and as in it's case, we can extend the 
linear region to cut the 'X' axis and look for the limiting initial angular displacement which 
does not show sharp exponential fall in oscillation time period. This works out to be 11 — 12°. 
Beyond this limit, the damping coefficient is large and a pronounced exponential fall is seen 
in the oscillation's time period (fig EI)- Kleppner and Kolenkow [7] have discussed the nature 
of (3 and have stated that it depends on the shape of the mass and the medium through 
which the mass moves. The amount of frictional force depends on the instantaneous velocity 
(d@/dt) of the pendulum {(3 being the proportionality constant). However, this nature of the 
frictional force (F = — /3dO/dt) is restricted for motion where velocity is not large enough to 
cause turbulence. Beyond angles of displacement of 11 — 12°, the frictional drag might not be 
following the linear relationship with velocity. This however needs further investigation. 

Before summarizing the results of our experiment, it would be of use to understand how 
the pendulum experiment is done in the undergraduate lab. The student records the time 
taken to complete 40-50 oscillation 8J oscillations from which the time period is calculated by 
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(degrees) 

Figure 9: Plot of T(/3 = 0) and T av vs angle. The solid line is the curve fit of T(f3 = 0) to 
eqnfJZJ) while the dash line is for visual aid to show variation of T av . 



dividing the total time taken by the number of oscillations measured in that time. We call 
this as T av (average). This term obviously does not take into consideration the influence of 
damping which is pronounced in large angle oscillations. This obviously leads to erroneous 
results. Figure (JHJ) compares T av and T(/3 = 0), the time period after accounting for damping 
with the pendulum's displacement. The figure clearly depicts the increasing disparity with 
large displacements. 

Thus, if a student performs the pendulum experiment without taking the necessary pre- 
caution of small angular displacement to get simple harmonic motion and in turn T Q , he or she 
would have to filter out the large angle correction and the damping coefficient. If no correction 
is made and 'g' is calculated using T av (listed in Table II), the variation in 'g' with angle 9, 
would increase (fig El hi)- The fall in time period with successive oscillations is evident in 
this experiment, since a micro-processor measures the time period. This would not have been 
evident in ordinary circumstances. Thus, the experimenter would not have been obvious of this 
and would only inco-operate large angle corrections, with no corrections regarding damping. 
The resulting variation is seen as the falling value of 'g' with angle in fig()10l ii). The true 
constant nature of 'g' (figElii) is only obtained when both corrections are inco-operated. 

Table II: Listed is the average time period T av that a student would measure manually 
along with the time period if he bothers to correct it for large angle oscillations. Also listed 
are the values of acceleration due to gravity he would have got with his time periods. 
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Figure 10: Plot shows how the calculated acceleration due to gravity, 'g', varies depending on 
what corrections have been done and it's variation with initial angular displacement. 



S.No. 





T 

- 1 - av 


Sav 


T a v 

i+e 2 /i6 


Scor 




(degrees^ 


(sec) 


(m/sec 2 ) 


(sec) 


(m/sec 2 ) 


1. 


5 


2.00477 


9.970 


2.00382 


9.979 


2. 


10 


2.00866 


9.931 


2.00486 


9.969 


3. 


15 


2.01246 


9.893 


2.00389 


9.978 


4. 


20 


2.00980 


9.920 


1.99464 


10.071 


5. 


25 


2.01563 


9.862 


1.99197 


10.098 


6. 


30 


2.02093 


9.811 


1.98696 


10.149 



Conclusion 

By doing the pendulum experiment with large angle displacements, calculations become com- 
plicated. As much as two informations have to be filtered out, the effect of large angle displace- 
ment and the damping factor. The damping coefficient is related to the initial displacement 
itself. These informations can only be processed if the time period of each oscillation is mea- 
sured. This is quite impossible manually and only a micro-processor interface is capable of 
highlighting these features. 
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Appendix 



The microprocessor program required for measuring the time period of eighty oscillations is 
listed below. 



Address 


Instruction 


Hex-code 


COOO 


MVI E 


IE 


C001 


160 D 


AO 


C002 


LXI H 


21 


C003 


00 


00 


C004 


CI 


CI 


C005 


MVI A 


3E 


C006 


00 


00 


C007 


OUT 


D3 


C008 


08 


08 


C009 


IN 


DB 


COOA 


09 


09 


COOB 


ANI 


E6 


cooc 


01 


01 


COOD 


CPI 


FE 


COOE 


00 


00 


COOF 


JZ 


CA 


C010 


15 


15 


con 


CO 


CO 


C012 


JMP 


C3 


C013 


09 


09 


C014 


CO 


CO 


C015 


IN 


DB 


C016 


09 


09 


C017 


ANI 


E6 


C018 


01 


01 


C019 


CPI 


FE 


C01A 


01 


01 


C01B 


JZ 


CA 


C01C 


21 


21 


com 


CO 


CO 


C01E 


JMP 


C3 


COIF 


15 


15 



15 



continued .. 



Address 


Instruction 


Hex-code 


C09n 


CO 


cn 


C!D21 


TXT R 


ni 

U X 


(T)99 

V_^UZjZj 




nn 

UU 


cn9s 


nn 


nn 

UU 




TNX R 


n^ 


Cn9"i 


TN 

XI i 


RR 


VJUZjU 


nQ 

u u 


nQ 

u u 




ANT 


Ffi 

XL/U 


cn98 


ni 

U J. 


m 

U X 




(HPT 


FF 

X XI/ 


V_ V Zji V 


ni 

u X 


m 

U X 


V. v ' z. i-f 


TZ 


CA 


Cn?C 

V U Zj V. 


94 


94 


Cn9R 

V_ V Z. _L_9 


nn 


cn 


Uaj -Li/ 


TNX R 


n^ 

uo 


Cn9F 


TN 

XI i 


RR 


cn^n 


nn 

Uc/ 


nn 

Uc/ 


firm 

VjU«J X 


ANT 


Ffi 

XL/U 


cn^9 


ni 

U J. 


ni 

U X 


cn^n 


CPT 


FF 

X XL/ 




nn 

UU 


nn 

UU 


cn^ 


T7 

J ZJ 


CA 


cn36 


2F, 

Zj J_J 


2F 

Zj J_J 


v V ■ ) I 


cn 


cn 


cn^s 


MOV M R 

ivxvy v xj> 


7n 

I u 


cn^n 


TNX H 

UMiV XX 


9^ 

ZjO 


cn^A 


MOV M C 


71 

1 X 


Q13R 


TNX H 


23 

ZjO 


cn^c 

V_y UO V_ 


RCR F 


1 R 

X 1.9 


C03D 


MOV A, E 


7R 


C03E 


CPI 


FE 


C03F 


00 


00 


C040 


JNZ 


C2 


C041 


21 


21 


C042 


CO 


CO 


C043 


HLT 


76 
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